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ABSTRACT 



A high frequency, broad-band transducer which has 
a broad radiation pattern in a half space is required 
for acoustic tracking of underwater vehicles. Two 
models of a flexural disk transducer, were built and 
tested.lt was found that the radiation loading in 
water completely damps the resonances due to standing 
flexural waves and the radiation appears to be due 
entirely to flexural waves propagating in the 
disk. Broad-band radiation patterns and broad bandwidth 
(60-80 kHz) were obtained. The discrepancies between 
theory and experiment are attributed in part to 
violation of the thin plate assumptions and to lack of 
knowledge of the actual velocity distributions. This 
approach offers promise for achieving the desired 
broad beam widt hs . However , additional development work 
is needed to obtain a transducer which has the curved 
face necessary for smooth hydrodynamic flow. 



4 



raSLE OF CONTENTS 



I. INTRODOCTION 9 

II. THEORETICAL CONSIDERATIONS 13 

III. EXPERIMENTAL PROCEDORES 19 

IV. SUMMARY OF AND DISCUSSION OF RESULTS 

AND SOME CONCLUSIONS 58 

7. RECOMMENDATIONS 62 

Appendix A-.COMPUTER PROGRAM FOR (0,2) MODE 63 

Appendix B:COMPUTER PROGRAM FOR PROGRESSIVE WAVE.. 65 

LIST OF REFERENCES 67 

INITIAL DISTRIBUTION LIST 68 



5 



LIST OF FIGURES 



1. Practical radiation pattern 12 

2. Calculated radiation pattern for normal mode 24 

3. Calculated radiation pattern for progressive wave... 25 

4. Transducer dimensions 26 

5. Transducer and driving element 27 

6. Admittance in air 28 

7. Displacement measurements 29 

8. Relative displacement amplitudes 30 

9. Amplitude distribution 31 

10. Amplitude distribution 32 

11. Admittance measurements 33 

12. Apparatus for beam pattern measurements 34 

13. Radiation pattern 35 

14. Radiation pattern 36 

15. Radiation pattern 37 

16. Radiation pattern 38 

17. Radiation pattern 39 

18. Radiation pattern 40 

19. Radiation pattern 41 

20. Radiation pattern 42 



6 



21. Admittance measurement 43 

22. Admittance measurement 44 

23. Admittance measurement 45 

24. Radiation pattern 46 

25. Radiation pattern 47 

26. Radiation pattern 48 

27. Radiation pattern 49 

28. Radiation pattern 50 

29. Radiation pattern 51 

30. Near field pressure amplitudes 52 

31. Near field pressure phase angles 53 

32. Near field pressure amplitudes 54 

33. Near field pressure phase angles 55 

34. Near field pressure amplitudes 56 

35. Near field pressure phase angles 57 



7 



ACKNOWLEDGEMENTS 



I am fortunate to have this opportunity to thank the 
people who have helped in the preparation of this 
thesis. First of all I am most indebted to my 
advisor, Prof essor O.B. Wilson, for his guidance and 
motivation. Also, I wish to express my sincere appreciation to 
Lieutenant Commander,?. Auns and Lieutenant A. Ford for their 
assistance in preparation of this thesis and Mr. S. Moeller 
for his assistance in constructing the transducers. 

I wish to acknowledge the unfailing support and 
tolerance of my wife,Bahar. 

I want to thank the Turkish and Dnited States Navies 
for creating this generous opportunity. 

Last tut not least I should also like to thank the 
Naval Torpedo Station Keypor t, Washington for their technical 
and financial support. 



8 




sot;- 



I 



INTR ODUCT ION 



For observing the position of submerged targets on 
underwater ranges, target-mounted sound sources are commonly 
used with fixed, bottom-mounted arrays. In order to 

insonify a large volume of water, essentially 
non-directicnal sources are required. 

Considering the geometry of one tracking range, that at 
the Naval Torpedo Station (NTS) ,Keyport, Washington, which 
has a typical vertical distance of 200 m and the distance 
between bottom- mounted hydrophone arrays of about 2000 
meters, the maximum slant range, r, can be found to be 
1020 m. Using the following formula for transmission loss 
(TL) , the transmission loss at the range, r, can be 

calculated: 



TL = 20 Log r + a.r 

For a sound frequency of 75 kHz, the attenuation 
coefficient, a, for sea water is about 0.03 dB per meter. 
The transmission losses are 91 and 52 dB at ranges of 1020 
and 200 m,respecti vely; hence, directional properties of the 
target-mounted sound source must compensate for this 40 dB 
difference in TL. Consideration of the pitch and roll 
motions of the underwater vehicle and the effect of the 
reflection from the surface has led to the specification 
that the source radiate uniformly within a conical volume 
150 degrees wide, centered on the vertical axis of the 
transducer. 
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The practical shape of the beam pattern is shown in 
Fig 1 has been calculated considering the difference in 
transmission loss with angle. 

In order to obtain the necessary source level at 75 
degrees from the main axis, it is required to drive the 
transducer at large amplitudes. This however, apparently has 
shortened transducer life (From informal communication with 
NTS) ; accordingly redistribution of total energy away from 
main axis was a major design consideration. It is felt that 
this would extend the Mean Time Between Failure (MTBF) of the 
transducer and associated amplifiers. 

Sufficient bandwidth for pulse code modulated telemetry 
and economy of the construction were additional goals. 

Formerly, this requirement was met by a pop-out 
transducer, utilizing a spherical piezoelectric ceramic 
element. However, for current high-speed targets the 
performance of the protruding radiator is unsatisfactory , due 
principally to fluid flow problems. Therefore, a number of 
designs of flush mount transducers have been built and used 
for this purpose. 

One approach to obtaining a broad beam radiation pattern 
is the use of a two element concentric piston transducer 
with amplitude and phase shading has been studied and which 
appears to offer pxomise (Ref. 4) .Another approach is that 
of Mr. M. Barlow, at the Naval Torpedo Station, who proposed 
the idea of and carried out some initial experiments using 
flexural waves in a small disk for this purpose ( private 
communication) . 

The intent of this study has been to explore further the 
use of transverse (flexural) vibrations in a plate as an 
approach in the development of a flush-mount transducer 
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capable of generating broad 
freguencies. 



beam patterns at high 
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Figure 1 - PRACTICAL RADIATION PATTERN 
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THEOSETICAL CONSIDERATIONS 



The derivation of the transverse wave equation for thin 
plates is given in texts on theory of elasticity (Ref . 1 p. 
214) . The equation is. 



where w is the displacement in the transverse direction, P 
is the density of the material, s is the Poisson’s ratio, Q 
is the modulus of elasticity (Young’ s Modulus) and h is the 
half thickness of the plate.lt has been assumed that h is 
small compared to the transverse wave length. 

For simple harmonic motion of a plate. 



, -rSn-ft 
w = Y (r , 0 ) e 

and a differential equation for Y can be written: 

Cv^- /) Y ^ 0 

where 

V i2nVVd-s>^) 



The solution of the first equation in polar coordinates, 
which is finite at r = 0,is, 
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Y = (m<p) 

S»n ^ 

where m is an integer. The solution of the second equation 
is obtained from the first equation by changing 
into i^ .It is called the Hyperbolic Bessel Function and 
is defined as 



I 

m 



(Z) 



.->n 

1 



•n 



(iz) 



Possible solutions for the simple harmonic oscillations 
of a plate are therefore given by the following expression: 



Y(r,^6) [a ^ 3 

A and B are derived from boundary conditions. 

The boundary conditions for a circular plate of radius 
a, clamped its edges are: 



The 



Y(a,0) = 0 and — = 0 at r=a. 

or 

first condition is satisfied by making 



B 



j JmCVal 
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and the second condition is satisfied by requiring ^ to 
satisfy 



I^Ofa)^ J^(j"r)-J^(j'a) ^ I^(^r) =0 at r=a. 



6f 



m"* ' m 



dr 



m 



Because If = Jf(f),this equation fixes the allowed 
frequencies. 



One can label the solutions for this second condition: 



where 






000 



A more complete tabulation of values of is given in 

Ref . 1, p. 215. 



The allowed frequencies for the normal modes are; 



f = 

'‘''I 2 



a 



.'/a 



seci-s^) 






2 . 2 



and the characteristic functions are: 






-r / TrrPmo ^ T f wT Frtry^ 

^ / Utt^) "V ^ i 
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On the other hand, the transverse wave speed for 
plate without losses is given at Ref 2,p. 508, 




where 



Q (Shf 
12 0-S^) 



is the stiffness of the plate for bending and M is the mass 
of the plate per unit surface. 



The following expression gives the value of the plate 
half thickness: 



h= 







a 




2. 3 



Hence, for a given frequency and plate material one can 
find a direct relationship between the thickness of the 
plate and the transverse wave speed c^j from equation 2. 3. One 
can also find a relationship between the radius a and 
thickness 2h by using equation 2. 2. In combining these two 
results the relationship between radius and thickness can be 
obtained. 

For aluminum plates (Ref 3) 

^ = 2700 kg / m^ 
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Q=7. 1 10 nt / 

s=.33 

For the (0,2) mode, and a frequency of 75 KHz 

i =0.2627 h 

cf = 1.47 10^ h 

b 

Computation of the radiation pattern function was based 
on a variation of a program worked out by Shaw (Ref 4) .In 
this case, the radiation pattern is the sum of radiations 
from a set of 20 concentric ring piston elements. 4n 
assumption is made about the amplitude and phase of the 
velocity of motion as a function of piston radius. 

The pattern for any ring, say number n, which has outer 
and inner radii of a^^ and a„_, respectively and a velocity 
amplitude U^, is computed by subtracting the 
ordinary , uniform piston radiation pattern function of radius 
a^_, with velocity amplitude Up. from that of a pistcn of 
radius a,^ and velocity 0^. 

A computer program shown in Appendix A was used to 
estimate the radiation patterns for combinations of various 
modes and radii. 

None of the patterns so obtained were found to resemble 
the specifications described before in the introduction 
section, but only (0,2) mode offered some promise. One of 
the calculated radiation pattern is shown in Fig 2. 

Also , radiation from circular progressive wave was 

estimated using the computer program given in Appendix B. 
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The equation for the displacement of an outgoing 
circular progressive flexural wave is: 



w=A af (k^r) 



where =Hankel function of the second kind 



Using the same method described above radiation patterns 
of various thickness and radius combinations were 
obtained. One of them is shown in Fig 3. 
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EXPE2IMENTAL PROCEDURES 



Based on the theoretical beam patterns, a plate with half 
thickness, h=1.5 mm and a radius r=2.33 cm, was chosen for 
a first try and the transducer shown in Pig 4 was 

constructed. Some dimensions were specified by the 

positioning and mounting hole for the transducer on the 
extender section of the test vehicle. 

To drive this transducer, a hollow piezo-ceramic 
cylinder, with radial polarization which has following 
specifications was used : 

Ceramic: GLENNITE HS21 

Length; 1.0 in 

Outer diameter: 0.5 in 

Wall thickness : 0.125 in 

Longitudinal resonance: 89.63 kHz 

Using a composite transducer configuration with mass-loading 
on the ends due to two pieces of aluminum (Pig 5), the 
longitudinal resonance was reduced to 73.6 kHz. In order 
to attach the masses to the element a nut and bolt ware 
used;however , for driving the plate the nut was found to be 
unnecessary . Without it, it was assumed that the longitudinal 
resonance would approach 75 kHz. 

Upon completion of the construction, electrical 
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admittance measurements of the transducer were made using a 
Dranetz Impedance- Admittance meter. The result for the 
unloaded transducer (in air) is shown in Fig 6. Primary 
resonances were found at frequencies: 

66.34 kHz, 73.26 kHz, 80.4 kHz, 105.2 kHz. 

The configuration shown in (Fig 7) and a Capacitive 
Displacement Probe (Ref 4 Fig 35) which averages 
displacement over the probe face were used to determine the 
modes of vibration for each resonance. After careful (and 
difficult ) measurements the amplitude distribution of 66.34 
kHz was found to resemble closely that expected for the 
(0,2) mode. The distribution at 73.26 kHz resembled somewhat 
that expected for the (0,3) mode. However the pattern was 
not symmetrical.lt is believed that this is one of the 
mixed-modes. Relative amplitude distributions are shown in 
Figs 8-10. 



Next the transducer was placed in the extender 
section. In this instance, admittance measurements showed 
small deviations from earlier resonance frequencies; 
however, when the mounted transducer was placed into 
water, none of the earlier resonances were found, except for 
the one at 72.27 kHz. Both admittance curves appear in Fig 
11. It appears that the radiation loading completely damps 
the normal modes of flexural vibration of the plate and that 
the resonance at 72.27 kHz is associated with the 
longitudinal vibration of the ceramic driver element. 

Measurements of the radiation pattern of the transducer 
were made in one of the Naval Postgraduate School* s anechoic 
water tanks using apparatus described below. The transducer 
was mounted in the wall of an aluminum cylinder which was 
the extender section of a 12-3/4 inches diameter torpedo. 



20 



The transducer baffle assembly can be rotated about a 
vertical axis, A potentiometer mounted on a shaft in the 
drive train provides an electrical signal which is 
proportional to angle of rotation which can be used to drive 
one axis of a graphical plotter. 

In order to eliminate the effects of surface 

ref lections , and minimize the effects of standing waves in 
the tank the pulse mode of operation was used. In this 
mode, the Boxcar Integrator, gated to the direct path pulse, 
generates a d.c voltage proportional to the amplitude of the 
pulse. This voltage is fed to a calibrated logarithmic 
converter which then generates the ordinate of the beam 
pattern in decibels. To improve signal to noise ratio (SNR), a 
high driving voltage (30 vrms or more) was used. The 

transducer and hydrophone were located with a two meter 
separation. The following test equipment was connected as 
shown in Fig 12: 

Hewlett Packard 7035 B X-y Recorder 
Tektronix RM45A Dual Trace Oscilloscope 
Hewlett Packard 465A Voltage Amplifier 
PAR 160 Boxcar Integrator 
Krohn-Hite DCA 50R Power Amplifier 
Harrison 6205 B D.C. Power Supply 
40 dB. Gate 

General Radio 1396A Tone Burst Generator 
General Radio 1192B Frequency Counter 
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Hewlett Packard 204C Oscillator 



Hewlett Packard 407A Power Amplifier 

LC-10 Hydrophone 

Envelope Detector 

Hewlett Packard 7561A Logarithmic Converter 

The radiation pattern for various frequencies are 
presented in Figs 13-21. 

It is believed that the measured radiation patterns are 
due principally to the progressive transverse wave in the 
plate and not, as first theorized, to the existance of 
standing flexural waves in the plate. 

The results show that radiation patterns below 60 kHz 
do approach the pattern which is shown earlier in 

Fig 1. However above 60 kHz, the level along the axis 
increases relative to that at larger angles. The broad beam 
width characteristics appear to change little with 
frequency. 

A second transducer was built with a thicker plate (4 mm) 
which should result in a larger transverse wave speed. 

Driving this new model (MOD 2) with the same ceramic 
element, the admittance curves shown in Figs 21-23 were 
obtained. The results were almost the same as earlier 
model. In water, the resonance was found at 72.24 kHz. 

Using the same configuration of test instruments shown 
in Fig 12, reasonably good radiation patterns were obtained 
below 78 kHz for this new model (Figs 24-29). 
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In an attempt to understand better the details of motion 
of the vibrating plate when in water, near field pressure 
measurements were made using the following equipment; 

Hewlett Packard 204C Oscillator 

General Eadio Frequency Counter 

Hewlett Packard 467A Power Amplifier 

Krohn Hite DCA 50R Power Amplifier 

Tektronix RM45A Oscilloscope 

LC-5 Hydrophone 

Dranetz Phasemeter 305 

Hewlett Packard 403B Voltmeter 

The LC-5 hydrophone was placed about 1 or 2 mm from the 
face of the transducer in the anechoic tank. A rack and 
pinion assembly supporting the hydrophone permitted moving 
the hydrophone across the face of the transducer for near 
field measurement £. Continuous excitation was used. The 
figures obtained for near field pressure amplitudes and 
phases are presented in Fig 30 and 31 for the MOD 1 
transducer, and in Figs 32 and 33 for the MOD 2 transducer 
when driven at the resonance frequency. Another measurement 
made at 60 kHz driving frequency gave the results shown in 
Figs 34-35. 
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Figure 3 - CALCULATED RADIATION PATTERN USING PROGRESSIVE 
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Figure 4 - TRANSDUCER DIMENSIONS 
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Figure 5 - TRANSDaCER AND DRIVING ELEMENT 
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Figure 6 - ADMITTANCE IN AIR 
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Figure 7 - DISPLACEMENT MEASUREMENTS 
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Figure 9 - AMPLITUDE DISTRIBUTION ON THE FACE OF THE 
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Figure 10 - AMPLITUDE DISTRIBUTION ON THE FACE OF THE 

TRANDUCER 
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Figure 12 - APPARATUS FOR BEAM PATTERN MEASUREMENTS 
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Figure 13 - RADIATION PATTERN 
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Figure 14 - RADIATION PATTERN 
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Figure 15 - RADIATION PATTERN 
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Figure 16 - RADIATION PATTERN 
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Figure 17 - RADIATION PATTERN 
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Figure 18 - RADIATION PATTERN 
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Figure 19 - HADIATION PATTERN 
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Figure 20 - RADIATION PATTERN 
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Figure 21 - ADMITTANCE MEASUREMENT 
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Figure 22 - ADMITTANCE MEASUREMENT 
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ADMITTANCE MEASUREMENT 
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Figure 24 - RADIATION PATTERN 
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azimuth angle 

Figure 25 - RADIATION PATTERN 
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Figure 26 - RADIATION PATTERN 
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F igure 28 - RADIATION PATTERN 
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Azimuth angle 

Figure 29 - RADIATION PATTERN 
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Figure 30 - NEAR FIELD PRESSURE AMPLITUDES 
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Figure 31 - N£AB FIELD PRESSURE PHASE ANGLES 
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Figure 32 - NEAR FIELD PRESSURE AMPLITUDES 



75.27 12Hz 
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Figure 33 - NEAR FIELD PRESSURE PHASE ANGLES 
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Figure 34 - NEAR FIELD PRESSURE AMPLITUDES 
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Figure 35 - NEAR FIELD PRESSURE PHASE ANGLES 



IV. 



SUMMI QZ and di^dssion of r esults and some 

CONCLUSI ONS 



There are some discrepancies between predicted normal 
mode frequencies and those actually measured which are 
difficult to account for quantitatively. They are most likely 
due to violations of the thin plate assumption and perhaps 
due to coupling between the modes of the plate with those of 
the longitudinal driver vibrator. 

From the admittance curves obtained when in water , the 
only resonance frequencies appear to be those of the driving 
element. 

This forces the conclusion that the radiation resistance 
completely damps the normal mode flexural resonances in the 
plate. The radiation which does occur then must be due to 
flexural waves propagating radially outward from the driver 
element located at the center of the plate. 

This conclusion is also supported somewhat by the 
results of the near field measurement s. In Figs 30,32,34, 
it is seen that the pressure amplitude is essentially 
constant over a diameter about equal to that of the driver 
element. Thus , the central part of the plate is behaving as a 
uniformly driven piston. The change in phase with radius 
along the plate beyond this center region corresponds to a 
wave speed which is higher than that in water. The phase 
speed of the flexural waves should increase with frequency 
and with plate thickness according to theory. In the three 
sets of data taken, the increase in phase speed with 
thickness was observed. The effect with frequency was 
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inconsistent and could not be explained. 

For the reason just stated the theoretical radiation 
patterns calculated based on the assumption of a propagating 
flexural wave from a point source should not be 
realized, because the source is not a point. 

The radiation patterns observed for the two models, 
shown in Figs 13-20 and Figs 24-29 can be summarized in 
the accompanying table(Table I) .The values of the angles to 
the outermost -6 dB points are tabulated at various 
frequencies. Any dip in response along the axial direction 
is also noted. 

Thera is a consistent assymmetry in the patterns which 
disappeared when a pattern was taken with a plane 
baffle. This is due to the location of the radiator nearer to 
one end of the cylindrical extender section used as baffle 
for these measurements. 

The angular widths to the -6 dB points are large, 
varying between about 120 to 148 degrees at frequencies 
which cover a band width more than 6 kHz on either side of 
resonance. These compare favorably to the 135 degrees width 
specified by the Naval Torpedo Station (fief. NAVORD 2814382 
Fig 2, Modified 2-28-73^ Flush 3-D Transducer Output 
Specifications) . 

The change in source level with frequency is not more 
than 4 dB in the band width between 60-80 kHz. 

The desired dip in source level at angles near the axis 
is not found at all frequencies. In fact, near resonance, 
there appears to be a small peak near the axis. It is not 
possible to account for these quantitatively . The reduction 
found in axial source levels is probably not broad enough to 
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be significantly useful in reducing overall power 
requirements. 

The maximum source level observed near resonance was 
about 75 dB ref. 1 microbar at 1 meter. This was for a 
200 volt rms drive voltage applied , which was the maximum 
undistorted value available with the equipment used in these 
tests. From the data on the electrical admit tance, the 
conductance at resonance is 700 micromhos. Thus at this 
drive voltage, the input power is about 28 watts during the 
pulse. 

According to private communications from NTS, the desired 
93 dB source level is achieved with about 400 watts input 
power during the pulse.lt is seen that this transducer , if 
driven with 400 watts, and if it were linear in 
response , would have a maximum source level of about 
87 dB, about 6 dB less than that desired. 
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V. 



R ECOM ME NDA TIONS 



Because of the promising results mentioned above , 
further study is recommended. 

The electromechanical efficiency can probably be 
improved somewhat if piezoelectric ceramic elements more 
suitable for this application were chosen. For example, one 
should use ceramic elements for the vibrator which have 
parallel field coupling and appropriate composition. The 33 
coupling instead of the 31 coupling used here should 
increase the electromechanical coupling coefficient and give 
a larger operating band width. The ceramic composition 
should be one designed for high power use (The Glennite HS-21 
used here is more appropriate for hydrophone use) . 

Some problems which must also be addressed if this 
approach is to be of practical use; 

The obtaining of a curved face to fit the curvature of 
the torpedo. This could be achieved by a rubber window which 
has same f*c as water which is molded to the flat disk. 

The radiation face of the transducer could be made 
larger with a possible resulting increase the source level. 

Further experimental tests on the effect of plate 
thickness with rubber window in place are recommended. An 
empirical approach is necessary to determine optimum plate 
thickness to get the desired radiation pattern. 
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DIMENSION A ( 23 ) , U ( 21 ) , PT Z( 20 ) , DARG ( 20 ) , X ( 30) , Y ( 30 ) , UDI F ( 20 ) 
1 ,FJ(20) 

R=RADIUS(M. ) 



APPENDIX A 



COMPUTER PROGRAM FOR NORMAL MODE SOLUTION 



OJ 

CD 

O 

OJ 

a 



in 

• 

CsJ 

f— ) 

• UJ 

o 

ox 

mo 



X 

m.-j 

• •* 

,-^UJ 

Luo: 

-o 

Xto 

CNJCO 

••LU 

- Qd 
lUO. 
Q^- 

o - 

OOX 



tm 

o • 

OCM 
O r— I 

moj 

m^ X 



O M 



a: 










(X 






LOO 


LU 










LU 




h- 


UJ^ 


1 — < 














o 


X ^ 


•• 










•• »— # 




h" 


X- 












^O^ 




+ 


OJ 


f 










• LUO 






C»-J 


•» 










*.>v •* 






<»— o 


o 






CQ 




O — 




Xw 


• rs.JXX 


Q 




t-H 


n 








— ^UL 


c\ji— <f <r 






•* 


cQ 




CNJ •— ow 




O— • 


^o.- - 


OC\J 


vO 


f— * 


< 




t-H-W- oo< 




1 Q 


^ •• •* 


•> •> 


*• 


II 


-+ 


O 


* 41 .-.^ 




— O 


1 XX 


o 


C\J^ 


X 


LL 


C\J 


^ CD-JV O 


o 




oorsjm 


II 


• II 


X 


(\JO 


CD •• 


*^oo 


CM 




j-o •• ^ 


ox 


0 *J<' *-< 




•»-00 


O — ^ 


O'-'OOO-^ 


«• 


o*-* 






(\J— 


V, 


LUCOLiJO II 


CsJ<C\JCO 1 MU 




II ^ 


luoo 


JLU f 


^ • 


CQOOl— • 


CQO<OCNJ.-h 


V, II V.OJOvOO 


II 





:cNjm 



H-CQOr^ 

^ II 

c\J— irocsj^:— I-J cos^rvi 

II osj II II o:a.-j o II 
UU^^N^O II <OcQ3: II II 0CQ0<a30C3rcQ 3300000— II <--CiOOO 
a3 0a.a.U.00Q<Xrv)XQ<0UL<Ul>^< <OU-OLU0Xa.cQO0 3:OF-O 



• >O<<Xc0'«’ <X II X— 

II II II '-H> II oj«— o II Lu<< 

II H- II H- II XCDJ-XH'LLr-l*-^-^'— -»• h-h- II nj + a. II II h-<H-XX 
OH-zu- II II II -irNio^— XH- *0.1— Xfsj«-<a:Qt: 

~ IIQXOOf-OlOO 

XZXF- 00 . 3 U- 11 . 



o 

o 



m 



C\J 



oo 

oo 

mm 



63 











cn 












— * t 












o 












^ f— ♦ 








iU 




M LL 








•— • 












• 


^- 


a. X 








r— 4 


o 


^ 00 








•o 


h- 


►— « •• 










< 


com 








-*-5C^ 


+ 


Q.Q • 








XC0< 


— o 


>-a-cNj 


o 








r 


cy^ •»»— 4 


*> 








w 


COi-^UJ 


o 




CO 




LL 


<1— •- 


m 






^Qw 


»— « 


WQ.X 


•» 




ir\ 


(/>Or>4 


Q 


O --CO 


> 




O' 


■«■ o:h- 


3 




p» 




CNJ 


— <Q- 


■S' 


O— ‘on 


X 


oo 


• 






OO*-* 




OrO 


f^ 




CNJX 


-iO - 


Q- 


ir\ •» 


m 


►<t c>OcQ 




<stx 


K 


•*ip^ 


V. 


— l-lf LU II U4f-40*-'U-l-»' •*'tXUJO 


O II 




• II 


rD II 0^1— 3 •^0'-Q3— 1 






a:^(X i£:xx<CLXO'--H-Q.2a. 


LU 


II 


II II w 


^ Q II 


*-*(NJLU< II 


•-H 




— •xi-mo-iot-'O II ^-^- II i-s:-' 


>- -J 



*-• CO*— • II O 

Qio II ^x^-oa:<<oo^-^-aQa:owo<^•z 
3Q-3XI— <Q<OQOQQ.<OCX3:tL>OOt>01^ 

ITi vO O 

o 



64 



APPENDIX B 



COHPaTER PROGRAM FOR PROGRESSIVE WAVS SOLUTION 
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